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A next-to-leading order QCD calculation of non-singlet spin structure function, 



^ • g± s is presented within the valon representation of Hadrons. In the valon model, it is 
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C*") . assumed that a nucleon is composed of three dressed valence quarks: the valons which 
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have their own internal structure, the valence quark with its associated sea quarks 
and gluons. The results are in good agreement with all available data from SMC, 
E143, HERMES and with the newly released data from COMPASS experiments. It 
appears that the small x tail of g± can be described by a single Regge-type exchange. 

X 

5h ' The relevant parameter of this exchange is given. Finally we show that the polarized 

proton structure function has a scaling behavior at small x. The relevant parameters 
of this behavior are given too. 
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I. INTRODUCTION 



Deep Inelastic Scattering (DIS) of leptons from the nucleon has served as an important 
tool for the investigation of the nucleon substructure and is one of the key areas for testing 
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the Quantum ChromoDynamics (QCD). 

Spin is a fundamental properties of the nucleon and the spin structure of nucleon has been 
the subject of heated debates over the past twenty years. The key question is that how the 
spin of the nucleon is distributed among its constituent partons. That is, the determination 
and understanding of the shape of quarks and gluon spin distribution functions, Sqj(x, Q 2 ), 
have became an important issue. 

We utilized the valon model jl| to study the polarized nucleon structure. In the valon model, 
it is assumed that a nucleon is composed of three dressed valence quarks: the valons. Each 
valon has its own internal structure, the valence quark with its associated sea quarks and 
gluons which can be probed at high enough Q 2 . At low Q 2 , a valon behaves as a valence 
quark. The valons play a role in scattering problems as the constituents do in bound state 
problems. It is assumed that the valons stand at a level in between hadrons and partons 
and that the valon distributions are independent of the probe or Q 2 . In this representation 
a valon is viewed as a cluster of its own partons. The evolution of the parton distributions 
in a hadron is effected through the evolution of the valon structure, as the higher resolution 
of a probe reveals the parton content of the valon. This model has yielded excellent results 
for un-polarized structure functions, l|-[6| . It has also been applied to the polarized nucleon 
structure function [7|, |8j with interesting results . 

In this paper we would like to concentrate on the non-singlet part of the polarized nucleon 
structure function; because of its simplicity and thus its transparency. In addition, there are 
more accurate data which are extended to fairly small x region: x < 0.01. That makes the 
comparison with the theory more meaningful. Recently COMPASS Collaboration released 
data on g^ s to test the Bjorken sum rule with more accuracy Therefore, in this paper we 
have further attempted to show the application of using this model for studying the nucleon 
structure functions. 

This paper is organized as follows. In section 2, a brief general outline is presented on the 
calculation of the polarized nucleon structure function in the valon model. Then we calculate 
the non-singlet spin structure function in section 3. Finally, section 4 is devoted to study 
the Regge behavior of g± s at small x and scaling behavior of g\. Then, we will finish with 
conclusions. 
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II. VALON MODEL AND POLARIZED HADRON STRUCTURE FUNCTION 

The connection between bound state problem at the hadronic scale that occurs only at low 
Q 2 and deep inelastic scattering at high Q 2 can be investigating by introducing the valons. 
Each valon is a dressed valence quark, i.e., each being a valence quark with its associated 
sea quarks and gluons which can be resolved only at high Q 2 . At low Q 2 valon behaves as 
a CQ because its internal structure can not be resolved. Thus the valon distribution in a 
hadron is the wave-function square of the CQs, whose structure functions are described by 
PQCD at high Q 2 Q. 




FIG. 1: (Color online) The schematic picture of the valon model. 

The valon model is essentially a two components model. In this framework, the structure 
function F 2 (x, Q 2 ) of a hadron is a convolution of the valon distribution G^ alon (y) and the 
structure function of the valon, F£ alon (z, Q 2 ) 

F 2 h (x : Q 2 ) = W 1 dyG h valon (y)Fr° n (-, Q 2 ) (1) 

valon ^ x y 

In a similar way the parton distribution functions in a hadron can be obtained by 

q(x, Q 2 ) = W 1 dyG h val0 Mr l ° n {j, Q 2 ) (2) 

valon Jx V 



where G^ alon (y) is the valon distribution inside the hadron. It means the probability of 
finding a valon with momentum fraction of y in the hadron. The description of the G^ alon (y) 
is given in jl, 2\. 
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valon(j) 






Pi 


aj 


h 


c i 


d, 


U 


3.44 


0.33 


3.58 


-2.47 


5.07 


-1.859 


2.780 


D 


-0.568 


-0.374 


4.142 


-2.844 


11.695 


-10.096 


14.47 



TABLE I: Numerical values of the parameters in Eq. (4). 



In polarized case, the helicity distributions of various partons in a hadron in the framework 
of this model are given by: 

• i 



5q*(x 1 Q 2 ) = J2 [ -^ valon {y)H al ° n 

j x y 



X 



-,Q 2 ) 
y 



(3) 



where SG^ alon (y) is the helicity distribution of the valon in the hosting hadron (probability 
of finding the polarized valon inside the polarized hadron). It is related to unpolarized valon 
distribution by: 



SG^y) = SFj^Gjiy) = N jy a *(l - yf'{\ + a jV ^ + b 3 y + Cj y Lb + d jV 2 ) 



(4) 



Gj(y) are the unpolarized valon distributions, where prefers to U and D type val- 
ons(Regrettably, the above equation was erroneous in Ref. 7J] in which 5Gj(y) was replaced 



by SFj(y) ). Polarized va 
for a number of hadrons 



on distributions are determined by a phenomenological argument 
7J, [8|. We summarized the parameters for Eq.(2.4) in Table 1. 




FIG. 2: Polarized valon distributions for U and D valons inside the proton. 
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We come back to Eq.(3), 5q^ alon (z = x/y,Q 2 ) is the polarized parton distribution in 
the valon . Polarized parton distributions inside the valon are evaluated according to the 
DGLAP evolution equation subject to physically sensible initial conditions. 



2tt ^/3o /v Nb± 2(3 ™ 



(5) 



5q s (n,Q 2 ) \ = { 3_ )s m» | a,(g 2 ) ^ L -(^)^p(°)" a,(Q§) £ _ ( ^ ) ^(o)n^ / <5 ? s (n, Qg) 



The detail of above equations are given at [10|. As presented in [3] and |8|, we have 
calculated the polarized nucleon structure function in the valon model. We have worked 
in MS scheme with Aqcd — 0.22 GeV and Qq = 0.283 GeV 2 . The initial motivation for 
this value of Qq comes from the phenomenological consideration that requires us to choose 
the initial input densities as S(z — 1) at Qq. The valon structure function has the property 
that it becomes 8(z — 1) as Q 2 is extrapolated to Qq (beyond the region of validity). This 
mathematical boundary condition means that the internal structure of the valon can not 
be resolved at Qq in the NLO approximation. It also means that at initial scale of Q^, 
the nucleon can be considered as a bound state of three valence quarks that carry all the 
momentum and the spin of the nucleon. As Q 2 is increased other partons can be resolved at 
the nucleon. It is also interesting to note that this value of Q 2 is very close to the transition 
region reported by the CLAS Collaboration. Measurement of the first moment of the proton 
structure function at CLAS shows that there is a transition region around Q 2 = 0.3GeV 2 



ll| . Therefore the initial input densities to solve the DGLAP equations inside the valon 



are 



5q NS (z,Q 2 ) = 5q s (z,Ql) = 6(z-l) (7) 

Sg(z,Q 2 ) = (8) 

Thus their moments are 
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Sq NS (n, Ql) = 5q s (n, Q 2 ) = [ z n 'H{z - l)dz = 1 (9) 

Jo 

5g(n,Q 2 ) = (10) 

In the valon model, the hadron structure is obtained by the convolution of valon struc- 
ture and its distribution inside the hadron. Having specified the various components that 
contribute to the spin of a valon, we now turn to the polarized hadron structure, which is 
obtained by a convolution integral as follows: 

^(*,q 2 ) = W 1 ^5G h valon (y)9r lon (^Q 2 )- (ii) 

valon Jx V y 

The valon structure is generated by perturbative dressing in QCD. In such processes with 
massless quarks, helicity is conserved and therefore, the hard gluons can not induce sea quark 
polarization perturbatively. According to this description, it turns out that sea polarization 
is consistent with zero [7[. This finding is supported by HERMES experiment and by 



released data from COMPASS experiment 12Hl5j. Therefore we have no sea polarization 
in our model. 

Using the initial conditions in Eq. (9) and Eq. (10), the calculation of the PPDFs 
inside the valon follow from the standard DGLAP evolution equations. The algorithm for 
calculation of the PPDFs inside the proton can be decomposed in the following three steps: 

• Calculating PPDFs in the valon by using the DGLAP equations; 

• With a phenomenological approach, one can find the helicity distributions for the 
valons, Figure 2. These functions are Q 2 -independent. Since we find the valon helicity 
distributions, one can use them to calculate the polarized nucleon structure up to 
Q 2 = 10 7 GeV 2 ; 

• By using the convolution integral (equations (3) and (11)), one finds PPDFs in the 
nucleon and polarized structure function; 



There is an excellent agreement between the model predictions with the experimental 
data for spin structure functions. A sample is given in Figure 3. 
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FIG. 3: (Color online) (a) Polarized proton structure function, xg\ at Q 2 = 5GeV 2 = 1-7V. (b) 
Polarized deuteron structure function, xqf at Q 2 = 3GeV 2 (jr = 1.33). The results from model 7\ are 



compared with the experimental data [9 
COMPASS. 
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mum- The data irom 



are newly released data from 



III. NON-SINGLET SPIN STRUCTURE FUNCTION 

The non-singlet polarized parton distribution function is defined as: Aq NS (x, t) = (Au 
Ad)(x,t); with t = ln^ . 
The DGLAP equation for Aq NS (x, t) is: 



d 
dt 



Aq NS = C d ^p(-)A q NS {y } t) 

Jx y \yj 



(12) 



where P(z = x/y) is the NLO spin-dependent splitting function 10|. As mentioned before, 
in our calculations Aq sea ~ , therefore Aq NS becomes (Auy — Ady)(x,t) which is shown 
in Figure 4 in comparison with other global fits. The non-singlet spin structure function is 
defined as: 



9i 



NS 



&-fi = 2\fi-g*/(l-1.5u D )] 



(13) 



where up = 0.058 accounts for the D-state admixture in the deuteron wave function. 

The g± s data is consistent with the quark model and the perturbative QCD predictions 
in the valence region x > 0.2 23]. 

In Figure 5(a) ,the x dependence of xg^ s is shown in comparison with data from HER- 
MES, E143, SMC and also with the newly released data from COMPASS [q]. The re- 
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x8uv-x8dv 




FIG. 4: (Color online) The non-singlet polarized parton distribution function in the valon model in com- 



parison with other global fits 



19|422 . 



suits are in very good agreement with experimental data for the entire measured range 
of x. In Figure 5(b) The evolution of the Bjorken integral, derived from the Figure 5(a), 
J^ min (^(xg^ s )dx = jx min {.9i S )dx as a function of x m i n is shown for the model compared 
with the recent HERMES and COMPASS Collaboration data [jj, Q. Note that about 50 
percent of the sum rule comes from x values below about 0.12 and that about 10-20 percent 
comes from values of x less than about 0.01. It shows that, g^ s receives a considerable 
contribution from the small x region. Thus, it seems that investigation of small x region of 
the structure function is important. In the following section we will consider g± s at this 
region. 

In Table 2 we compared the integral over different x ranges at different scales of Q 2 , as de- 
termined from the valon model with the experimental results from COMPASS, HERMES, 
.^143, E154, E155, SMC and JLAB also with the recent results from NN Collaborations 



24 



27] . We have used data from JLAB for Q 2 > O.bGeV 2 to make sure that 



the non perturbative effects are small 
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x 9i 



0.1 

0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
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• HERMES 2 6 
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© E143 
AAC 
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MODEL 




Integrate [glNS, {x, . 9} ] 
(b) 



0.2 



0.15 



n 



□ COMPASS200 9 

O HERMES200 6 

A MODEL 



A 

2 



A' 



online) (a) xg^ s at Q 2 = 5GeV 2 compared with the experimental data and the results 



FIG. 5: (Color 

from global fits [I9|-|2l| = 1.25) . (b) The integral of 6 over the range 0.02 < x < 0.9 measured by 
HERMES and COMPASS Collaborations as a function of the low x limit of integration, x m i n , evaluated at 
Q 2 = lOGeV 2 in comparison with our results. 

IV. REGGE BEHAVIOR OF gf s AND FULL g\ AT SMALL X 

In all the results from experimental data for unpolarized and polarized structure func- 
tions, it is seen that these structure functions increase when x decreases and Q 2 increases 
for fixed values of x and Q 2 respectively. The small x behavior of spin dependent structure 
functions reflects the high energy behavior of the polarized virtual compton scattering total 
cross section with increasing total CM energy squared W 2 since W 2 = Q 2 {\ — 1). When 
W 2 » Q 2 , x is small and W 2 ~ Q 2 /x and then at this region the structure functions 
have scaling behavior. This is, by definition, the Regge limit and so the Regge pole ex- 
change picture is therefor quite appropriate for the theoretical description of this behavior 
. The small x or high energy behavior of the spin structure function is an important 



issue for the extrapolation of data needed to test spin sum rules for the first moment of g\ . 
The small x measurements, besides reducing the error on the first moment, would provide 
valuable information about Regge and QCD dynamics at small x where the shape of g± is 
particularly sensitive to the different theoretical inputs. 

In the case of unpolarized structure function, F2, it is believed that a Regge trajectory com- 



bined with a soft and a hard pomeron can accurately represent the experimental data 



29] 
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Experiment 


x — Range 


o 2 


1 i 


This Analysis 


COMPASS 


0.004 < x < 0.7 


3 


0.175 ±0.009 ±0.015 


0.1421 


HERMES 


0.021 < x < 0.9 


5 


0.1479 ±0.0169 


0.1381 


E143 


0.03 < x < 0.8 


2 


0.149 ±0.016 


0.1276 


E143 


0.03 < x < 0.8 


3 


0.164 ±0.023 


0.1301 


E143 


0.03 < x < 0.8 


5 


0.141 ±0.013 


0.1327 


E154 


0.03 < x < 0.8 


5 


0.168 ±0.010 


0.1327 


E155 


0.03 < x < 0.8 


5 


0.176 ±0.008 


0.1327 


SMC 


< x < 1 


10 


0.198 ±0.023 


0.1626 


SMC 


< x < 1 


5 


0.174 ±0.024 ±0.012 


0.1569 


JLAB 


0.001 < x < 0.8 


0.592 


0.1027 ± 0.0228 ±0.0052 


0.123 


JLAB 


0.001 < x < 0.8 


0.707 


0.0945 ±0.0201 ±0.0151 


0.1298 


JLAB 


0.001 < x < 0.8 


0.844 


0.1021 ±0.0193 ±0.0174 


0.134 


JLAB 


0.001 < x < 0.8 


1.20 


0.1307 ±0.0192 ±0.0145 


0.1323 


JLAB 


0.001 < x < 0.8 


1.44 


0.1522 ±0.0186 ±0.0089 


0.1433 


JLAB 


0.001 < x < 0.8 


1.71 


0.1605 ±0.0182 ±0.0069 


0.1318 


JLAB 


0.001 < x < 0.8 


2.05 


0.1678 ±0.0177 ±0.0056 


0.1475 


JLAB 


0.001 < x < 0.8 


2.44 


0.1666 ±0.0167 ±0.0045 


0.1492 


JLAB 


0.001 < x < 0.8 


2.92 


0.1789 ±0.0106 ±0.0035 


0.1511 


NN Collaboration 


0.021 < x < 0.9 


5 


0.1315 ±0.0144 


0.1381 



TABLE II: Comparison of the integral over different x range at different scale of Q 2 , as determined 
from the valon model, with the experimental results from COMPASS, HERMES, E143, E154, 
E155, SMC and JLAB also with the results from NN Collaborations. 



It is interesting to investigate this issue in the polarized case. 

The Regge pole model gives the following small x behavior of the structure functions 



g\(x,Q 2 



28) 



g[(x,Q 2 )= li (Q 2 )x^ (14) 

where g[(x,Q 2 ) denote either singlet (gl(x, Q 2 ) = gl(x,Q 2 ) + g™(x, Q 2 )) or non-singlet 
(gi s (x, Q 2 ) = gl(x, Q 2 ) —g™(x, Q 2 )) combination of structure functions. Therefore we expect 
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that describe the small x behavior of g^ s and g\ structure functions with one and two Regge 
exponents respectively. It appears that the present g^ s data for available small x in the 
interval 0.0001 < x < 0.01 can be described with a single Regge type exchange as : 

g* s = g\- g\ ~ Ax aR ^ (15) 

The Regge intercept which governs the small x physics has smooth Q 2 dependence, in 
which case one would see g± s rising like x a where —0.5 < a < also at low Q 2 and in 



the measured "small x" region 30| . This value varies from -0.13 to -0.3 when Q 2 is moved 



from Q 2 = 2GeV 2 to Q 2 = lOGe^ 2 in the valon model. According to the results of Ref. 31], 
asymptotic scaling of g± s depends on one variable Q 2 /x 2 only, instead of two variables x 
and Q 2 with the constant intercept equal to 0.42: 

g? s ~ (Q 2 /x 2 ) Ans/2 ~ Q Ans x~ Ans ;A ns = 0.42 (16) 

However, it is valid for a very small x only. The applicability region of that analysis 
is x < 10~ 6 . Figure 6 shows the non-singlet spin structure function at Q 2 = 5GeV 2 
for small x region. In Figure 7, we fit the non-singlet spin structure function at small x 
(0.0001 < x < 0.01) and find the associated aR egge : 

g N.s = Ax <*r w *. A = 0.173, aRegge = _o.323 (17) 



A. Scaling behavior of cf x at small x 

Having analyzed the small x behavior of gi S , it is interesting to see how g\ behaves as 
x — > 0. Because the valon model has very good agreement with existent small x data, it is a 
good candidate to show the small x behavior of g\ at small x. Actually the results for this 
scaling behavior should be compared with small x data for g\. 

We show that the polarized proton structure function has this scaling behavior for 1.2 < 
Q 2 (GeV 2 ) < 100 at small x (10~ 5 < x < 10~ 2 ) 

2 

g p 1 (x,Q 2 ) = J2^(Q 2 )x £ * (18) 

i=l 

Where a« and Si are constants and the functions fi{Q 2 ) have this simple general form 



12 




FIG. 6: (Color online) Small x behavior of xg± s at Q 2 — 5GeV 2 calculated by using the valon model in 
comparison with the experimental data . 




lxlCT 4 2xl(T 4 5xl0- 4 0.001 0.002 0.005 



x 

FIG. 7: Small x behavior of g± s at Q 2 — 5GeV 2 calculated by using the valon model and using the best 
fit to calculate the Regge exponent. Data points are from the model and the goodness of fit is: x 2 = 0.99. 



HQ 2 ) = (j^Qi) Di (19) 
HQ 2 ) = HQ 2 )g(Q 2 ) (20) 

(21) 
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Where 

g(Q 2 ) =9o + 9iLog(Q 2 ) + g 2 Log(Q 2 ) 2 + g 3 Log(Q 2 ) 3 (22) 

The results for parametrization of g\ are summarized in table 3. (We should noted that 
if we try to do a fit with only one regge exponent and only with the first term in Eq. 18, we 
give a fit with x 2 — 0.99. So, It is clear that we can describe the polarized proton structure 
function with two exponents, x 2 = 0.996, better than one exponent ) 

0.008 
0.006 
0.004 
0.002 


-0.002 

10 s 10* 10" 3 10" z 

X 

FIG. 8: (Color online) Small x behavior of xg\ at Q 2 = 5GeV 2 calculated by using the valon model and 
using the best fit to calculate the scaling exponents at small x . The experimental data selected for a various 
of Q 2 at small x region. 




In Figure 8 we show the Small x behavior of xg\ at Q 2 = 5GeV 2 calculated by using the 
valon model and using the best fit to calculate the scaling exponents, £j. 
As a result, we concluded that only by using two scaling exponents, we can describe the small 
x behavior of g\ (x, Q 2 ) well. This behavior has also been seen in the small x tail of F%(x, Q 2 ) 



structure function 



28] . The question has been raised 32|, |33| whether the observed rise on 



g± (x, Q 2 ) follows from one or two pomeron exchange ( a polarized analogue of the one or 
two pomerons question). However, looking at ^(x,Q 2 ), it indeed requires two pomerons! 
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Parameters 


values 


£l 


0.196 


£2 


0.094 


ai 


0.0215 




0.0513 


A 


0.759 


£>2 


0.434 




17.538 


5i 


-11.809 


52 


2.652 


53 


-0.200 


Qo 


1.300 


X 2 (Goodness of fit) 


0.996 



TABLE III: global fit parameters obtained by fitting the Eq. (4.4) with the data extracted from 
the valon model at small x. 

V. CONCLUSION 

In this paper we calculated the non-singlet spin structure function, gi S , of the nucleon 
in the valon model. The results of these calculations are in excellent agreement with all 
experimental data for the entire measured range of x. We also study the small x behavior 
of non-singlet spin structure function and the Regge behavior of g± s to calculate the Regge 
exponent, ctR eg ge- Finally, we studied the scaling behavior of g\ at small x. We conclude that 
only by using two scaling exponents, we can describe the small x behavior of g\ well. This 
is very similar to existence of two soft and hard pomeron to describe the small x behavior 
of F 2 (x, Q 2 ) in unpolarized case. It is shown that the valon model can predict the polarized 
nucleon structure functions for the entire measured x range very well. The validity range 
for using the valon model is 0.5 < Q 2 (GeV 2 ) < 10 7 and 10~ 6 < x < 1. In this model sea 
quarks polarization is consistent with zero (the finding that was confirmed by very recent 
experiments at HERMES and COMPASS), so the polarized hadron structure functions can 
build only by polarized valence and gluon distribution by finding only one type of polarized 
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valon distribution for each kind of valons. Between 0.3 < Q 2 (GeV 2 ) < 0.5 we should consider 
other effects such as combination of resonance physics and vector-meson dominance that are 
important at low Q 2 . Chiral perturbation theory may describe the behavior of polarized 
nucleon structure close to threshold. These issues will be considered in our future works to 
improve our model. 
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